An important objective of evolutionary biology is to understand the processes that govern phenotypic variation in natural populations. We assessed patterns of morphological and genetic divergence among coastal and inland lake populations of nine-spined stickleback in northern Sweden. Coastal populations are either from the Baltic coast (n = 5) or from nearby coastal lakes (n = 3) that became isolated from the Baltic Sea (< 100 years before present, ybp). Inland populations are from freshwater lakes that became isolated from the Baltic approximately 10,000 ybp; either single species lakes without predators (n = 5), or lakes with a recent history of predation (n = 5) from stocking of salmonid predators (~50 ybp).

Coastal populations showed little variation in 11 morphological traits and had longer spines per unit of body length than inland populations. Inland populations were larger, on average, and showed greater morphological variation than coastal populations. A principal component analysis (PCA) across all populations revealed two major morphological axes related to spine length (PC1, 47.7% variation) and body size (PC2, 32.9% variation). Analysis of PCA scores showed marked similarity in coastal (Baltic coast and coastal lake) populations. PCA scores indicate that inland populations with predators have higher within-group variance in spine length and lower within-group variance in body size than inland populations without predators. Estimates of within-group P_ST (a proxy for Q_ST) from PCA scores are similar to estimates of F_ST for coastal lake populations but P_ST >F_ST for Baltic coast populations. P_ST >F_ST for PC1 and PC2 for inland predator and inland no predator populations, with the exception that P_ST <F_ST for body size in inland populations lacking predators.

Baltic coast and coastal lake populations show little morphological and genetic variation within and between groups suggesting that these populations experience similar ecological conditions and that time since isolation of coastal lakes has been insufficient to demonstrate divergent morphology in coastal lake populations. Inland populations, on the other hand, showed much greater morphological and genetic variation characteristic of long periods of isolation. Inland populations from lakes without predators generally have larger body size, and smaller spine length relative to body size, suggesting systematic reduction in spine length. In contrast, inland populations with predators exhibit a wider range of spine lengths relative to body size suggesting that this trait is responding to local predation pressure differently among these populations. Taken together the results suggest that predation plays a role in shaping morphological variation among isolated inland populations. However, we cannot rule out that a causal relationship between predation versus other genetic and environmental influences on phenotypic variation not measured in this study exists, and this warrants further investigation.